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An effective low-energy description for multi-leg spin-1/2 Heisenberg ladders with an odd number 
of legs is proposed. Using a newly developed Monte Carlo loop algorithm and exact diagonalization 
techniques, the uniform and staggered magnetic susceptibility and the entropy are calculated for 
ladders with 1, 3, and 5 legs. These systems show a low-temperature scaling behavior similar to 
spin-1/2 chains with longer ranged unfrustrated exchange interactions. The spinon velocity does 
not change as the number of legs increases, but the energy scale parameter decreases markedly. 
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Recently, antiferromagnetic Heisenberg spin-1/2 lad- 
der attracted much interest, following the discovery of a 
spin gap in the 2-leg-ladder [EJ . Also the crossover from 
the single chain to the two-dimensional (2D) square lat- 
tice, obtained by assembling chains to form "ladders" of 
increasing width, is far from smooth |2| . Heisenberg lad- 
ders with an even number of legs, nj, have a spin gap 
and short range correlations, while odd-leg ladders have 
no gap and power-law correlations. These theoretical 
predictions have been verified experimentally, in mate- 
rials such as (V0) 2 P 2 Oy H and the homologous series 
of cuprates Sr n ^\Cu n j r \02n PI, which contain weakly 
coupled arrays of ladders. 

Here, we concentrate on odd-leg ladders. Our goal is to 
derive a low-energy description in terms of S=l/2 chains 
with longer range effective interactions, and examine the 
evolution with increasing number of legs, m. The Heisen- 
berg Hamiltonian for ladders is 
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where i and j enumerate the rungs, r, r' label the legs, 
and the sum marked by <-> (\) runs over nearest neigh- 
bors along legs (rungs). Periodic boundary conditions 
are chosen along the leg direction and open boundary 
conditions perpendicular to it. For the known materi- 
als we expect the superexchange to be roughly isotropic 
(J_i_ = J||). However, it is educational first to consider 
the strongly anisotropic limit (Jj_ ^> Jn). 

In the completely anisotropic limit (Jii/J_l = 0), each 
eigenfunction is a direct product of one- rung states whose 
lowest-lying multiplet is a spin doublet, separated by a 
gap of order from the first excited state. The ground 
state of the whole system is therefore 2 L -fold degener- 
ate. A finite value of Jy lifts this degeneracy. Our goal 
is to formulate an effective Hamiltonian, H ff, in this 
2 L -dimensional subspace M. of rung doublets which de- 
scribes the low-energy properties. For the case of the 
3-leg-ladder, to third order in Jii/Jx, we get 
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where S/ ' = Sj.i + Sj, 

j th rung, J n = J_LEA a n,A(j^-) 
103/243, a 2 ,i=0, 



the total spin of the 

with aj. i=l, ai,2= - 
a 2 , 2 =-8/27, a 2)3 =-49/162, 

3/T2 



1/9, ai, 3 

a 3 ,i=a 3 , 2 =0, a 3 , 3 =32/243, and J = (16/81) Jf/J_ 

The last term in Eq. (||) will be neglected, since the 
corrections in the energy of the low- lying energy states, 

(3) 

due to this term are small. H^, ff has then the form 
of a single chain with effective nearest neighbor (n.n.) 
coupling Ji, next-nearest neighbor (n.n.n.) coupling J 2 , 
and with exchange coupling J 3 between rung spins sepa- 
rated by three unit cells. Therefore the low-lying energy 
states of the 3-leg-ladder can be mapped onto those of 
a Ji — J 2 — J 3 chain. In this effective system the n.n.n. 
interaction are F (J 2 < 0), while the n.n. n.n. interactions 
are AF (J 3 > 0). Consequently, both the second and the 
third term in Eq.(|^) enhance the overall AF quasi- long 
range order. Note, that the third order corrections in 
Eq.(||) affect 3 1 and J 2 strongly since the corresponding 
coefficients are large. So one must perform the calcula- 
tions at least up to third order. 

(3) 

To test Hg ff , we calculate the temperature dependent 
uniform susceptibility, x(T), for 3- leg ladders || using 
the Quantum Monte Carlo (QMC) loop algorithm 
and compare to susceptibilities obtained for the effective 
Ji — J 2 — J3 model. We consider large enough systems, 
such that finite size effects are negligible. All results are 
extrapolated to a Trotter time interval At — > 0. Further, 
we compare with recent results obtained by Greven et al. 
§ using the same algorithm. 

At low temperatures, where only the states in M. are 
relevant, the susceptibilities of the 3-leg ladders with 
small J||/Jj_ coincide with those of the corresponding 
Ji — J 2 — J 3 chain. This can be seen in the inset of Figju 
where we show the susceptibility per rung of the 3-leg 
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ladder with Jii/Jj. = 0.2 together with the susceptibili- 
ties of the corresponding effective models in first, second 
and third order in Jii / Jj_. While the first order effective 
model (a single chain with only n.n. interactions) gives 
only a qualitative description of the 3-leg ladder at low 
temperatures, x(T) of the effective model in third order 
in Ji|/J_i_ coincides with the susceptibility per rung of the 
3-leg ladder up to a crossover temperature. Above this 
temperature, x(T) of the 3-leg ladder is larger, due to the 
presence of additional states in the 3-leg ladder which are 
not included in the subspace A4. 
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FIG. 1. Susceptibility of the 3-leg ladder for different J || /Jx 
and of the corresponding effective models. The filled symbols 
show the data for the 3-leg ladders and the open symbols 
those of the corresponding third order effective model. The 
crosses show the susceptibility of the corresponding Ji — J2 
chains in the mapping of the 3-leg ladders to Ji — J2 chains 
(for details see text). The inset shows the susceptibility per 
rung of the 3-leg ladder with J||/J_l = 0.2 together with those 
of the corresponding effective models in first, second and third 
order in Jy/Jx- The error bars are smaller or of the order of 
the symbols. 

The crossover temperature is of the order of the gap 
A to the higher lying states in the 3-leg ladder. It can 
be estimated best by considering the entropy of both the 
3-leg ladder and the corresponding effective model (Fig. 
^]). Just above the crossover temperature, the additional 
states lead to a rise oc e -A / T in the entropy. Conse- 
quently, fitting the form e _A / T to the difference of the 
entropy per rung of the 3-leg ladder and the entropy of 
the corresponding effective model gives an rough esti- 
mate of A. We find A rj Jj_ more or less independent of 

JJ/J.L. 

From Figjlj, it is clear that with increasing J||/J_l the 
quality of the description of the 3-leg ladder by the third 
order effective model becomes worse. The qualitative 
features of the temperature dependence of the suscepti- 
bility, however, are correctly given even in the isotropic 
case J||/Jj_ = 1, e.g. the slope of x(T) is increasing, while 
the zero-temperature value x(0) remains more or less 



constant, as Ji|/J_i_ increases (see Fig.| 
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FIG. 2. Entropy of the 3-leg ladder for different J|i/Jj_ 
(solid lines) and of the corresponding effective models (dashed 
lines). 

Calculating the effective Hamiltonian of the 3-leg lad- 
ders to n th order leads to a spin-chain Hamiltonian with 
interaction terms between spins which are separated by 
up to n unit cells. All these interactions are invariant 
under translations and rotations. Consequently, for very 
low temperatures the 3-leg ladders can be mapped onto 
the k = 1 Wess-Zumino-Witten (WZW) nonlinear a- 
model |^,|l^|, which is determined by a spinon velocity, 
v, and a energy scale parameter, To. For T <C To, x(T) 
of this model reads @Q up to: O ((In T)" 3 ) 



X(T) = h ~t~ 
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ln(ln(T /T) + 1/2) 



ln(T /T) 
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X(T) approaches its T = value x(0) = l/(27rv) with 
infinite slope. To should be < A, and characterizes the 
interactions between the spinons. The smaller To, the 
stronger the interactions, and the faster x(T) increases 
with temperature. 

The low temperature regime of the universality class of 
spin chains with a rotationally and translationally invari- 
ant Hamiltonian (to which also the 3-leg ladders belong) 
is determined by only two parameters, v and To. How- 
ever, the determination of v and To for the 3-leg ladder 
is difficult. QMC calculations cannot be performed down 
to low enough temperatures such that a fit to the above 
form (||) gives reliable estimates for v and To . Especially 
To is considerably underestimated in all cases. 

To overcome this problems we first map the 3-leg 
ladder to a Ji — J2 chain and then study the one-to- 
one mapping of this Ji — J2 chain to the WZW model: 
(Ji, J2) <-> (v, To). The mapping to Ji — J2 chains is 
always possible since the low-T range is characterized by 
only two parameters which can be chosen as Ji, J2 in- 
stead of v, T . The mapping of the 3-leg ladder to Ji — J 2 
chains is done as follows. For small Jii/Jj. we can use the 
third order result for Ji and J2 (Eq. (0)), neglecting J 3 , 
since J3 is small. Otherwise we fit x(T) of Ji — J2 chains 
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for low T to x(T) of the 3-leg ladder (see Fig. [l]) which 
gives estimates of the values Ji and J2 (see Table ||). 



6.0 



•— • finite size scaling 

— 0(3) non linear o-model 




FIG. 3. Spinon velocity, v, and energy scale parameter, 
To, of the Ji — J 2 chains, determined by finite size scaling 
analysis (circles) respectively by mapping to the 0(3)-non 
linear a-model. 

The mapping (Ji, J2) <-> (v, Tq) can be studied, using 
exact diagonalisation methods, v and T are determined 
by the finite size scaling of the energy gap between the 
excited state E(k = tt, S z = 1) and the ground state 
E(k = 0, S z = 0) @: 
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tt, S z = 1)-E(k = 0, S z = 0) = 

ln(ln(L/L ) + 1/2) 
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where E(k, S z ) is the lowest energy with wavevector k 
and z-component of spin S z for a chain of length L. Ln is 
the characteristic scaling length of the chain. As a conse- 
quence of the equivalence of the imaginary time direction 
and the space direction, Ln = v/Tn. Fitting Eq.(Q) for 
different lengths using exact diagonalisation gives v as 
well as To. The results for different fractions J2/J1 are 
plotted in Fig.||. An alternate possibility to determine v 



and To is to map the Ji — J2 chains onto the 0(3) non- 
linear rj-model with 9 = tt. For the coupling constant 
g one finds g = 4 (1 - 4J 2 /Ji)~ 1/2 @. Therefore the 
spinon velocity can be written as: 
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Logarithmic corrections, on the other hand, depend on 
a mass scale A, which is generated dynamically in the 
tr-model , and for small g, A = 1/a exp(— 2n/g), where 



a is the lattice spacing |13|. Since Tq oc A, 



T n oc exp 



4J 2 
Ji 
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The results for v and T (Eq.(§) and Eq.(§)) are plot- 
ted in Fig|] together with those obtained by finite size 
scaling. 

In Table |, the estimated values v, To (following the 
above described procedure) are given for various 3-leg 
ladders. The spinon velocity increases, while To de- 
creases with increasing Jm/Jx- 

The values for v and To can now be put back into 
Eq.(||) which gives the low temperature susceptibilities 
of the corresponding 3-leg ladders. This is shown at the 
example of J M /3± = 0.2 in Fig.f§ Together with the QMC 
data for the 3-leg ladder itself and its third order effective 
model, we get the susceptibility for the 3-leg ladder with 
good precision on the whole temperature range from zero 
temperature to high T. 
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TABLE I. For low temperatures, 3-leg ladders can be 
mapped onto Ji — J2 chains. The corresponding coupling con- 
stants, Ji the 3-leg ladders are also listed. The value enclosed 
in bracket was obtained by fitting the £-data |9j to Eq. (]7J). 
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FIG. 4. Susceptibility of the 3-leg ladder with Jy/J± = 0.2 
and its effective models. The error bars are smaller or of order 
of the symbols. 

Greven et al. recently calculated the correlation length 
£ of isotropic ladders || , using also the QMC loop algo- 
rithm. For T <C To the inverse of the correlation length 
in the WZW-model can be written as H] 

1 ~ T ( 2 1 . Hn(ln(T /T) + l/2) \ 

e(T) ~ V ln(T /T) + 2 ln 2 (T /T) J 1 1 
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Fitting the data |J for ni = 3 to the above form (Eq.(0)) 
also gives an estimate of To. In the isotropic case we find 
To = 0.34 Ji| which is in good agreement with our result 
(see Table @). For the 5-leg-ladder T is already < O.Uy. 

Finally, we examine the static structure factor C(n, tt), 
defined by (N s = number of sites) 



C(k x ,k y 



- E 



a ik x (i-j)+ik y (r-r') 



< Sj T Sj T i > 



and show that it can be calculated from the effective 
model without introducing additional physical parame- 
ters. C(7r,7r) can be written as £ y (-l) i+j < SfSf > 

/3L, where Sf = S^i — S;^ + Si, 3 is the staggered spin of 
one rung. In the limit J|| = the correlations of the stag- 
gered rung-spins is simply related to those of the uniform 
rung-spins by: 



< b; b: >i^j= A < b; bj >i^j, 



(8) 



with A = 25/9 for all temperatures, where only the states 
in M. are relevant. For J|| ^ the "wave function" of the 
spin-1/2 degree of freedom at each rung is spread out over 
a certain number of unit cells. On one hand, this gives 
rise to the longer range interactions, described above, and 
on the other hand, it affects the correlations < SfSr > 



and < S[ ot Sj ot 



>. For large distances |i— j| ^s> 1, however, 

< Sf Sf > and < S] ot S* ot > are still related by Eq.(§) 
but with a renormalized value of A. 

Here, we calculate A numerically, using QMC simula- 
tions. The values are plotted in Fig. || for different Jy/J± 
and extrapolating to J|| = leads to a value A = 2.78 
which agrees with the analytical result 25/9. At low 
T, the correlation length £ ^S> 1 and C(7r,7r), as well as 
C(tt,0) = EijC- 1 ) 1 ^ < S? ot S* ot > /3L are dominated 
by correlations with |i — j| 3> 1. Additionally, C(ir, 0) 
per rung of the 3-leg ladder differs only slightly from 
C(tt) of the single chain. Consequently at low temper- 
ature, C(7r, 7r) 3 _ log = AC(7r, 0) 3 _ic g ~ AC(7r) singlo ch ./3. 
Greven et al. have calculated G(ir, 7r)3_i eg for J|| = 
and C(7r) S i ng i c h. 0. From their data we determine 
C(7r, 7r) 3 _i cg /C(7r) s i„ g ie ch. = 2.55, which is in good agree- 
ment with our value A/3 = 2.64 (see Fig|j). 

The above considerations can be generalized to an ar- 
bitrary odd-leg ladder with ni legs. There are no qualita- 
tive differences. The overall AF quasi-long range order, 
however, increases with increasing in. Therefore, espe- 
cially the ratio IJ2/J1I of the effective Hamiltonian H c ff 
is larger. The logarithmic corrections increase markedly 
(To decreases) and the T — > behavior sets in at lower 
temperature as ni increases. The zero temperature value 
x(0), on the other hand, is almost independent of nj. 
This implies that the spinon velocity of odd-leg ladders 
depends only slightly on the number of legs. 

In conclusion, we have proposed an effective spin-1/2 
chain model which describes the low-energy properties 
of odd-leg ladders. The temperature dependence of x, 



£, and C(7r,7r) for the effective system is shown to be 
consistent with that of the original model at low enough 
temperatures. The effective model requires two param- 
eters, e.g. a spinon velocity v and an energy scale, To. 
With an increasing number of legs, v does not change, 
but To decreases rapidly. The exchange interactions of 
the corresponding effective model become longer ranged, 
and antiferromagnetic correlations are enhanced. 

The calculations were performed on the Intel Paragon 
of the ETH Zurich. The support of the Schweizerischer 
Nationalfond is gratefully acknowledged. 




FIG. 5. The fraction A =< Sf Sf > / < Sf ot S- ot > at 



|i — j| 3> 1 for different Jy/Jj 



[1 
[2 
[3. 
[4 

[5 

[6 

[7 

[8] 
[9] 

[10; 

[11 
[12 



[13 
[14 



E. Dagotto, J. Riera, and D.J. Scalapino, Phys. Rev. B 
45, 5744 (1992). 

For a review, see E. Dagotto and T.M. Rice, Science 271, 
618 (1996). 

D.C. Johnston, J.W. Johnson, D.P. Goshorn, and A.J. 
Jacobsen, Phys. Rev. B 35, 219 (1987). 
Z. Hiroi, M. Azuma, M. Takano, and Y. Bando, J. Solid 
State Chem 95, 230 (1991); M. Azuma et al., Phys. Rev. 
Lett. 73, 3463 (1994). 

B. Fris chmuth, B. Amnion , M. Troyer, Phys. Rev. B (in 
press), 3ond-mat/9601025. 

H.G Evertz, G. Lana, and M. Marcu, Phys. Rev. Lett. 
70, 875 (1993). 

H. G. Evertz, M. Marcu, Nucl. Phys. B (Proc. Suppl.) 30, 
277 (1993). 

U.J. Wiese, and HP. Ying, Z. Phys. B 93, 147 (1994). 
M. Greven, R.J. B irgeneau, and U.-J. Wiese, preprint, 
cond-mat/9605064 

S. Eggert, I. Affleck, and M. Takahashi, Phys. Rev. Lett. 
73, 332 (1994). 

K. Nomura, Phys. Rev. B 48, 16814 (1993). 

I. Affleck, Les Houches Lecture Notes, in: Field, Strings, 
and Critical Phenomena, ed. E. Brezin and J. Zinn- Justin 
(North Holland, Amsterdam, 1988). 

G. Sierra, Jour. Math. Phys. A29, 3299 (1996). 

K. Nomura and M. Yamada, Phys. Rev. B 43, 8217 

(1991). 



4 



